Mott-insulator phases of spin-3/2 fermions 
in the presence of quadratic Zeeman coupling 
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We study the influence of the quadratic Zeeman effect in the Mott-insulator phases of hard- 
core spin-3/2 fermions. We show that contrary to spinor bosons, any quadratic Zeeman cou- 
pling preserves a SU(2) ® SU(2) symmetry, leading for large-enough quadratic Zeeman coupling 
to an isotropic pseudo-spin- 1/2 Heisenberg antiferromagnet. Depending on the scattering lengths, 
on ID lattices the quadratic Zeeman coupling can induce either a Kosterlitz-Thouless transition 
between a gapped dimerized spin-3/2 phase and a gapless pseudo-spin- 1/2 antiferromagnet, or a 
commensurate- incommensurate transition from a gapless spin-liquid into the pseudo-spin- 1/2 an- 
tiferromagnet. Similar arguments allow to foresee corresponding transitions on ladder type and 
square lattices. We analyze various observables which should reveal in experiments these phases. 
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Ultra-cold gases in optical lattices constitute an ex- 
traordinary tool for the analysis of strongly-correlated 
gases under extremely well-controlled conditions [l], 0] , as 
highlighted by the observation of the superfiuid to Mott- 
insulator (MI) transition in ultra-cold bosons Q , recently 
followed by the realization of the metal to MI transi- 
tion in two-component fermions 0, HJ. Due to super- 
exchange, the MI phase of spin-1/2 fermions is expected 
to acquire a magnetic Neel (antiferromagnetic) ordering, 
whose observation is the goal of active on-going efforts [6| . 

Optical traps permit the simultaneous trapping of 
various Zeeman sublevels, allowing for multi-component 
(spinor) gases. Spinor bosons have attracted a large in- 
terest due to their rich ground-state physics and spinor 
dynamics Q. Recently, spinor fermions are also at- 
tracting a rapidly-growing attention, motivated by ex- 
periments on BEC-BCS crossover in two-component 
fermions [2[ and the availability of multi-component 
fermions, including three-component Li gases [8], four- 
component 40 K [9(, spin-3/2 fermions as 135 Ba and 
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Ba [ll|], and Fermi-degenerate Yb [K 



Multi- 
component fermions present a wealth of novel phases. 
Pseudo-spin- 1 fermions allow for color superflui dity 
and trions, whereas spin-3 / 2 gases are even richer [l3 - 16|. 
including quartets. Ultra-cold Yb opens the fascinating 
possibility of 5 , C/(6)-symmetric spin-5/2 gases [l7|. 

A rich physics is also expected for repulsive spin-3/2 
fermions [l4l . Il5 |. which in quarter filling may undergo 
a MI transition. Whereas for spin-1/2 a Neel phase 
is expected, the MI in spin-3/2 presents a richer mag- 
netic structure given by a gapless spin-liquid or a gapped 
dimerized phase, depending on the interatomic interac- 
tions. Contrary to the spin-1/2 case, where spin-changing 
collisions are absent and the quadratic Zeeman effect 
(QZE) is irrelevant, the QZE is crucial for higher-spins, 
as shown in spinor condensates Q. However, in spite of 
its experimental relevance, the QZE has been typically 
ignored in the analysis of the magnetic properties. 
In this Letter, we explore, for the first time to our 
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FIG. 1: Phase diagram for spin-3/2 fermions with QZE. We 
depict the chirality r (density plot), the singlet-triplet cross- 
ing in the excitation spectrum (squares), the jump of the ex- 
ponent 7 to —1 (circles), the critical q c expected from two- 
band theory (triangle), and the J2/J = 0.25 line resulting 
from a strong-coupling analysis (white curve). See text. 



knowledge, the role of the QZE in the magnetism of 
spin-3/2 MI phases. We show that contrary to spinor 
bosons [Hj], any QZE preserves a SU(2) ® 5(7(2) sym- 
metry which results at large QZE in a pseudo-spin- 1/2 
isotropic Heisenberg antiferromagnet (iHAF). Depending 
on the scattering lengths, our results suggest for growing 
QZE a Kosterlitz-Thouless (KT) phase transition (PT) 
from a gapped dimerized spin-3/2 phase into the iHAF, 
or a commensurate-incommensurate (C-IC) PT between 
a gapless spin-liquid again into the iHAF. We analyze 
in detail observables which may characterize the phases 
(Fig. [TJ and PTs in future experiments [19( ■ 

We consider a balanced mixture of spin-3/2 fermions 
in a ID lattice [20] with N m = N- m (N m = number of 
fermions with spin projection m), such that the magne- 
tization M = J2 m mN m = 0. For a deep lattice and low 
filling (single-band regime) the system Hamiltonian is 



H = -t 
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at site 
with 



annihilate fermions with 

{F] Ttx f \ttl\ . Tfi^) the Clebsch-Go: 
n m ,j — V'mjV'mj- t is the hopping rate between adja- 
cent sites, and 50,2 = ^irh 2 aQ y2 /M characterize the s-wave 
channels with total spin and 2 (the only available due 
to symmetry), with ao t2 the scattering lengths and M 
the atomic mass. Although typically ao and a 2 are simi- 
lar, their values may be variated by means of microwave 
dressing (2fj or optical Feshbach resonances [23 |. Below 
we use G = (flrj+fl^A and g = (g 2 - 5o)/(ff2 + go)- Tne 
interactions preserve A4 and hence the linear Zeeman ef- 
fect induced by a magnetic field does not play any role 
for a fixed M.. However spin-changing collisions (which 
re-distribute the populations of the different components 
while preserving Ai) are crucial. Finally, the QZE, char- 
acterized by the externally controllable (by means of a 
magnetic field or microwave or optical dressing) constant 
q, is crucial in spinor gases, and plays a key role below. 

For large-enough interactions, G>f, we may consider 
the hard-core case, with maximally one fermion per site. 
For a chemical potential fi larger than a critical /Lt c (t) the 
system enters into the MI regime with one fermion per 
site. We focus below in the magnetic properties of these 
MI phases. Within the MI regime we perform perturba- 
tion theory obtaining an effective Hubbard Hamiltonian 
with super-exchange interactions, which for two adjacent 
sites {1,2} is of the form H\ 2 = Hq 



H sc , where 
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contains self-energies and effective hoppings without 
spin-changing, whereas a simultaneous hopping and spin- 
changing is given by H sc = c sc (b\ /2 B 3/2 + H.c.J , with 

B m = i>m,%i>-m,\ - i>-m,ui>-m,i- We introduce above: 
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with r)± = ±1, X± = G — 5q/2 ± [V + g 2 G 2 ] 1/2 , and 
tan</> = ([V +g 2 G 2 } 1 / 2 + 2q) /gG. 

Due to the spin-changing collisions, the single-site 
two-particle eigenstates with zero magnetization are not 
states with spins m and — m (|m, — m)), but the linear 



combinations |+) = cos <p\ - 3/2,3/2) + sin cj>\ - 1/2, 1/2) 
and |-> = -sin^|-3/2,3/2)+cos^|-l/2, 1/2) (with re- 
spective eigenenergies X± ). This mixing leads to H sc and 
the non-trivial q and g dependence of c 3 / 2 , Ci/%, and c sc . 
Note that C3/2, C1/2, c sc may diverge at specific values 
for which spin-changing collisions and QZE enter in res- 
onance. Below we consider only non-resonant situations 
for which the hard-core formalism remains valid. 

For q = the grou nd-state of H = J2i 0Iu y 
depends on g 1J, [l5j. For — 1 < g < the ground- 
state is a gapless spin liquid phase (with 3 gapless spin- 
modes). This phase includes the exactly solvable 5C/(4) 
line (g = 0, i.e. g = g 2 ) [23]. For < g < 1 the ground- 
state is a dimerized (spin-Peierls) phase which exhibits a 
spin gap and finite long-range dimer-dimer correlations 
lim»^oo(AA+n) = fo, with Di = {-l)%(S i - 1 + S i+1 ), 
with S the spin-3/2 operators. At g = the system 
undergoes a KT-like transition between both phases. 

The ±1/2 and ±3/2 manifolds are just linked by H sc . 
Hence, at large-enough q > 0, when e = c sc /{ci/ 2 — 
C3/2 — 2q) <C 1, the ±3/2 manifold is favored, and the 
system acquires a pseudo-spin-1/2 character. Below we 
restrict to q > 0, but similar reasonings apply to q < 0, 
for which the ±1/2 manifold is favored. Projecting onto 
the ±3/2 manifold and up to second order in e > 0, 
the Hamiltonian reduces to an isotropic Heisenberg form, 
H12 ~ JSi • S2, where J = — 2(c 3 / 2 — 2ec sc ) and S de- 
notes pseudo-spin-1/2 operators S z = {h 3 / 2 ~ n_ 3 / 2 )/2, 
5+ = ^ , 3 / / 2 V-'-3/2 an( l S- = 5J|_. Since J > for all 
<70!<?2) 1 > 0, the system reduces to an iHAF for large- 
enough q. This is true at any order in e due to a hidden 
SU(2) symmetry (see below). This must be compared 
to the case of spinor bosons where virtual transi- 
tions between manifolds lead for large q to an anisotropic 
Heisenberg Hamiltonian, demanding a fine tuning of the 
microscopic constants to map the system to an iHAF. 

Interestingly, spin-3/2 fermions, which show a hid- 
den 50(5) symmetry for q — [13j . retain at any q a 
SU(2) (g> SU(2) symmetry, generated by a direct prod- 
uct of two 5/7(2) spin algebras operating in the ±1/2 
and ±3/2 manifolds, respectively. The SU{2) genera- 
tors are 5^ = 1/2 JZii^i - n~ a ,i), S% = J2i/ l P a ,i'>P-<x,i, 
with a = {1/2,3/2}. They belong to the 50(5) alge- 
bra generating the 50(5) symmetry for q — and also 
commute separately with the QZE term. Note that this 
hidden SU(2)<E> SU(2) symmetry is related neither to the 
hard-core constraint nor to quarter filling, being rather a 
generic feature of 50(5) fermions with QZE. This sym- 
metry might be very helpful for future numerical simula- 
tions on four-component fermions in magnetic fields. For 
large q this symmetry results in the mentioned iHAF. 

Growing q induces PTs between the q — phases and 
the iHAF. We have obtained the ground state for vari- 
ous q and g by using the Density-Matrix Renormalization 
Group (DMRG) method of Ref. [H, with 36 sites, open 
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boundary conditions, and matrix dimension 20. The 
asymmetric phase diagram (Fig. [T]), discussed in details 
below, is characterized by a re-entrant dimerized phase. 

Case of g > 0. Bosonization shows that the dimerized 
phase is robust at small q, which induces a finite chirality 
t= i[{N 3/2 + N_ 3/2 )- (N 1/2 + N_ 1/2 )} Kq = qG/t 2 (for 
q <C 1). At large q bosonization is not appropriate, and 
one must instead descend from the iHAF decreasing q. 
Spin-changing processes lead to an AF frustrating next- 
nearest neighbour exchange J 2 (~ q~ 2 ) between pseudo- 
spins-1/2, resembling a frustrated spin-1/2 Ji — J 2 AF 
chain, which presents a PT at J 2 /J\ ~ 0.25 (25|. Hence a 
PT can be anticipated between the iHAF and the dimer- 
ized phase since when lowering q and increasing g the 
ratio J 2 /J increases. Fig.Q]shows the curve J 2 /J = 0.25 
obtained from perturbation theory, which is in good 
agreement with the results discussed below. 

A strong-coupling analysis is however not a reliable 
proof of the existence of an actual PT, and numerical 
calculations must be performed to confirm it. In a finite 
chain the two ground states of the dimerized phase (de- 
generate in the thermodynamic limit) split into a unique 
ground state and an excited one separated by an energy 
gap exponentially small in the system size. Thus, for fi- 
nite chains the lowest excited state in the dimerized phase 
is unique (for < q < q cr ). In contrast, the lowest ex- 
cited state above the Heisenberg ground state (q > q cr ) 
for a finite size chain is a degenerate triplet. If the PT 
is of KT type, as in a frustrated spin-1/2 J\ — J 2 AF 
chain, a level crossing between the lowest excited singlet 
and triplet states should occur. Exact Lanczos diagonal- 
ization results, for up to 12 sites with periodic boundary 
conditions and performing finite size scaling, confirm in- 
deed this crossing. The extrapolated q cr lies in the region 
expected from the DMRG calculation of r (grey region 
on the g > side of Fig. [1]). Although this certainly 
suggests a direct KT PT from the dimerized phase to 
the iHAF, for large g > 0, our size limitations are too 
severe for small g -C 1, and there we cannot exclude the 
existence of intermediate phases (26[. Hence, the iHAF 
chain undergoes a dimerization transition. For decreas- 
ing q the dimerized phase of the pseudo-spin- 1/2 chain 
adiabatically connects with the dimerized phase of spin- 
3/2 fermions, as they share similar properties. 

Although beyond the scope of this Letter, a rich 
physics is also expected in ladder-like and square lattices, 
since 0(q~ 2 ) terms induce frustrating diagonal (J x > 0) 
and next-nearest neighbor ( J 2 — 0.5J X ) exchanges. On 
ladders an Ising PT is expected from a rung-singlet phase 
to a columnar-dimer state with growing frustration [30j, 
while on a square lattice the Neel state may undergo a PT 
into a columnar-dimer state 31] possibly via a second- 
order PT showing deconfined quantum criticality (32l |. 

We are interested in the ground-state characterization 
at finite q for g ^ 0. The dimer phase is characterized 
by a finite fo and a singlet-triplet gap A in the magnetic 
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FIG. 2: Chirality r (a) and exponent 7 (b) as a function of q, 
for different values of g. DMRG results with 36 sites. 



excitation spectrum [27j . Both fo and A vanish for the 
iHAF. We monitor as well spin-spin correlations, (S l z ■ 



SI) oc \j — i| 7 [281 ], and the chirality t, which ranges 
from (q = 0) to 1 (pure iHAF). Fig. da) shows DMRG 
results for r. For weak g > the system presents still a 
quasi-saturation of r, which smoothly converges to 1 for 
large g. 7 behaves similarly, jumping abruptly to —1 for 
small g > (Fig.[T]). fo and A evolve similarly. 

Case of g < 0. In this case there is a PT between 
two gapless phases for increasing q. For g = (no spin- 
changing collisions) the critical properties as a function of 
q resemble those of a two-band model [2^] , where atoms 
with to = ±1/2 (±3/2) act as fermions at the lowest 
(second) band, and the QZE difference 2q resembles the 
band gap. For q < qo = 2 In 2 ~ 1.386 the magnetic order 
is suppressed due to "orbital" effects, and the system has 
three massless spinons. On the contrary for q > qo the 
"orbital" degeneracy is lifted, and the manifolds ±1/2 
and ±3/2 completely decouple in the ground state. When 
this occurs, r saturates to 1, and the system reduces, as 
mentioned above, to a pseudo-spin-1/2 iHAF with J = 
2t 2 /G. Therefore at qo there is PT from a gapless spin- 
liquid into a gapless AF pseudo-spin-1/2 chain. At the 
phase transition the exponent 7 jumps to —1 (29j . We 
have determined at g = the critical value q ~ 1.35, in 
good agreement with the expected value. Our results are 
consistent with 1 — r ~ Vqo~~q when approaching the 
PT for growing q. Hence at g = and qo there is a C-IG 
PT 33] between the two gapless phases. Moreover, g = 
and qo is a multicritical point for three phases: spin- 
liquid, pseudo-spin-1/2 iHAF, and the dimerized phase. 

The region g < smoothly connects with g — since 
perturbations from g = into g < are (margin ally ) ir- 
relevant in the renormalization group (RG) sense 15( and 
the symmetry dynamically enlarges to SU(A). One could 
hence expect that the g < region behaves similarly to 
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the g = case for growing q. There is, however, an im- 
portant distinction, since for g 7^ r is not a good quan- 
tum number, never saturating for finite q. The situation 
resembles that of electrons in an external magnetic field 
in the presence of spin-non-conserving processes, which 
are irrelevant in the effective description but explicitly 
break the microscopic symmetries of the lattice model 
spoiling magnetization as a good quantum number [jyj. 
A plausible scenario for g < is that the QZE induces a 
C-IC PT, so that chirality-non-conserving processes re- 
main irrelevant all the way, and do not modify the nature 
of the PT which takes place at g — 0. Numerical simu- 
lations must be performed to confirm this scenario. The 
PT cannot be identified by studying r, and instead we 
study 7 [iilj]. A jump in 7, if present, will confirm the 
C-IC nature of the PT. At small g < we obtain that 
the PT retains the main features of that at g = 0, and 
that indeed 7 abruptly jumps into —I (Figs. [T]and[2£b)). 

In summary, the QZE strongly modifies the MI phases 
of hard-core spin-3/2 fermions. Interestingly, contrary to 
bosons, a large SU (2) <S> SU{2) symmetry remains at any 
q, resulting at large q in a gapless pseudo-spin- 1/2 iHAF. 
For g > our results suggest a KT PT from a gapped 
dimerized spin-3/2 phase into the iHAF. On the contrary, 
for g < there is a PT between a gapless spin-liquid into 
the iHAF, which (at least for low \g\) belongs to the C- 
IC universality class. The phase diagram (Fig. [T]) has 
a non-trivial asymmetry characterized by a re-entrant 
dimerized phase. We have studied various observables, 
as chirality and spin correlations, which may reveal these 
phases in future experiments. Note that spin-changing 
collisions provide a unique opportunity to study exper- 
imentally a field-induced C-IC PT in four-component 
fermions, contrary to two-component ones where a mag- 
netic field cannot induce a C-IC PT due to the conserva- 
tion of magnetization. Note also that a richer physics is 
expected for higher spins, e.g. 5/2, where we expect for 
growing q a sequential transition from a spin-5/2 into a 
pseudo-spin-3/2 and finally a pseudo-spin- 1/2. 
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